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Abstract 

We report the analysis of radial characteristics of the flow of granular 
material through a conical hopper. The discharge is simulated for various 
orifice sizes and hopper opening angles. Velocity profiles are measured 
along two radial lines from the hopper cone vertex: along the main axis of 
the cone and along its wall. An approximate power law dependence on the 
distance from the orifice is observed for both profiles, although differences 
between them can be noted. In order to quantify these differences, we 
propose a Local Mass Flow Index that is a promising tool in the direction 
of a more reliable classification of the flow regimes in hoppers. 


1 Introduction 


The study of granular flows has attracted great interest from physicists and 
engineers due to its applicability and rich phenomenology [l]-[^, particularly 
due to the applications in the industrial and mining processes where it plays a 
crucial role [^. Moreover, in the last decades, with the increasing computational 
performance of codes 

based on the Discrete Element Method (DEM), there has been a new surge 
of interest on the subject with the possibility of studying large systems at the 
scale of real experiments. One of the oldest and most important cases of granu¬ 
lar flow is the discharge of grains through a bottleneck. Efforts have been made 
to understand the multitude of features presented by this type of flow such as 
clogging [5}|^, density fluctuations [S , segregation [^[^, strong fluctuations 


on the forces acting on silo walls 11 


12 and the collapse of silos during dis¬ 


charge 13 . Another known property of granular discharges is that two distinct 
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flow regimes can be observed — mass flow and funnel flow —, depending on 
the material properties of the grains and walls, and mostly on the hopper ge¬ 
ometry (grain-wall friction, opening angle, orifice size, among others) 17 


In engineering applications, mass flow is usually desired, since it is the flow 
where grains leave the hopper in a “flrst-in, first-out” order. Oppositely, funnel 
flow is characterized by the presence of stagnation zones, so that part of the 
grains is slower or immobile and acts as a funnel next to the hopper walls. In 
some applications, such as typically the agricultural ones, the stagnant material 
may deteriorate causing direct loss of product as well as efficiency. Therefore, 
the better understanding of these phenomena, particularly the influence of the 
design of silos, can impact directly the industry of bulk solids. 

The flow regime classification is often based on visualization of the steady 
state. Through the decades it has evolved from the sole visual observation 
of stagnation zones on experiments to the time averaged velocity spatial 
distributions, either by color scaled two-dimensional fields or profile curves along 
a certain linear path. Velocity profiles, particularly, are widely applied in the 
field of granular discharge flow 1^ - 21 . The funnel flow regime is then recognized 


by the presence of stagnation zones in the velocity fields or a large difference 
between the velocity in the central region and that close to the walls. In order 
to quantify this difference, an index is defined 14 the Mass Flow Index 
corresponding to the ratio of the average particle velocity next to the 
hopper wall, to the average particle velocity at the hopper centerline, 
that is, Jmf = It is defined such that if Jmf <0.3 the regime is funnel 

flow, and mass flow otherwise. 

Ketterhagen et ai used this index to create design charts for wedge 
shaped and conical hoppers, that is, they measured Jmf from simulations for 
different pairs of grain-wall friction and hopper wall angle and plotted them in 
a phase diagram to serve as a guideline for the designing of hoppers. For this, 
they averaged and Vc over a large region spanning along the hopper walls 
and the centerline, respectively. 

The velocity dependence on the vertical coordinate is usually assumed to be 

Tips 


equal to that of an incompressible fluid 


In the case of an incompressible 


fluid flowing through a conical hopper with an angle a with the horizontal plane, 
we can easily obtain the velocity dependence on the vertical coordinate. Being 
the volume flow rate defined as Q = Av, where A is the cross sectional area and 
V is the velocity. The flow rate must be constant along the direction of the flow, 
provided that the density is homogeneous. Thus, consider a point at height 
on the hopper opening, where Qo = and an arbitrary point at height 2 ; in 

the hopper, where Q = Av, then 


R. 

Qo 


Av TT {z tan a)‘^ V 


^o'^o IT {zo tan a)‘^ Vo 


= 1 


and, from that, 


V = { — 


-2 


( 1 ) 


Therefore, for an incompressible fluid, the velocity decays as a power law along 
z with exponent 2. However, if Equati on holds for a granular discharge, the 
procedure used by Ketterhagen et al. \2‘^ may lead to an Jmf value that is 
highly sensitive on position and size of the averaging region. This would cause 
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their design charts to be applicable only if the same specifications were used. In 
fact, Gonzalez-Montellano et al 18 calculated the ratio of the wall and center 


velocities, but averaging them over small regions at a few different heights along 
a square cross-section silo with a hopper-bottom that also has square cross- 
section. Although the ratio is approximately equal to one all along the vertical 
silo, they observed that inside the hopper the ratio largely fluctuates reaching 
values below 0.3 at some points, even though they consider the discharge to be 
in the mass flow regime. 

Granular flow may differ from incompressible fluids as there are voids and, 
therefore, density heterogeneities occur, even though each grain is roughly in¬ 
compressible by itself. In fact, one can expect that a larger density would 
be observed at positions farther from the hopper outlet, where the velocity is 
smaller and the grains are in a densely packed structure, as compared to re¬ 
gions closer to the orifice, where the velocity is larger and grains are looser. 
As shown by Sielamowicz and Czech 24 in experiments discharging amaranth 


seeds through a two dimensional hopper, the assumption of an incompressible 
fluid is not always valid. They found discrepancies from the 2D equivalent to 
Equation specially, close to the outlet. Additionally, their flow regime was 
closer to funnel flow as there is the formation of a channel of flowing grains in 
the central region of the hopper, while the grains close to the walls of the hopper 
are stagnant. 

In order to contribute to this discussion, we present in this paper the velocity 
profiles along the centerline and along the walls of a conical hopper and we show 
that they behave roughly as power laws. Nevertheless, there are some differences 
between Vc and profiles that carry information about the flow regime. We 
propose the analysis of the Local Mass Flow Index profile as a method capable 
of extracting that information in a clear way. Also, this method should be 
applicable independently of the validity of any assumption on the velocity radial 
dependence. Furthermore, we believe that the Local Mass Flow Index analysis 
may lead to a more precise definition of the flow regime boundaries. 


2 Model 


The microscopic model used for the grain interactions was described in details 
by Poschel and Schwager 25 . A particle i, with radius Ri and at position r^. 


is in contact with a particle j if ^ = [{Ri Rj) — \rj — ^|] >0, where ^ is the 
overlap length. The normal force fn is given by the Hertzian Law: 


fn = V ^eff 




( 2 ) 


where ^ is the normal relative velocity; Feff = {E- ^ Ej ^ is the effective 
elastic coefficient; Ei = (4/3)[Ti/(l — where Yi is the Young modulus and 
Ui is the Poisson ratio of the material of particle i; i^eff = 

effective radius and A is a normal damping constant [^. We also impose pos¬ 
itive values for so that unphysical attractive forces are avoided. Therefore, 
the normal force vector is 


Fn = max(0,/n)en, 


(3) 


where = {vj - fl) / \fj - rl|. 
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The tangential force is implemented in the spirit of Cundall and Strack’s 
seminal paper 27 . In the static regime the force is represented by an elastic 


term, analogously to a spring being attached to the points of the first contact 
of each of the particles surface, and a dissipative term. 


fo = -ktC-ltVt (4) 

where ( is the spring displacement, kt is the spring stiffness, 7 ^ is a tangential 
damping constant and Vt = Vij — en{en • Vij) is the tangential velocity. Only 
the tangential displacements are considered, so that ( and Vt are parallel to 
each other. If the tangential force is larger than the static friction limit, that 
is /o > jJisfni then the particles are in the dynamic regime, where the force is 
proportional to the normal force. 


fc — k'dfni 


(5) 


being fis and respectively, the coefficient of static and of dynamic friction. 
Therefore, the tangential force is 


f /o if /o < l^sfn (static regime) 

1 fc if /o > l^sfn (dynamic regime) 


We use material parameters as follows: E = 10^ Pa, corresponding to T ^ 
10^ Pa and u ^ 10“^, = 0.8, fid = 0.6, kt = 100 N/m, 7 ^ = 0.1 Kg/s, A = 
0.01 s. In order to avoid crystallization, the radii of the spheres are randomly 
chosen from a uniform distribution from an interval of 5% deviation around 
R = 5 mm. The density is 8 x 10^ kg/m^, therefore each grain weights about 
4.2g. Additionally to the contact forces, gravity acts on each particle and the 
acceleration of gravity in the simulations is set to ^ = 10 m/s^ in the direction 
of —z. The integration of Newton’s equations of motion was performed using 
a fourth order Gear’s predictor-corrector algorithm [^. The time step was 
chosen as At = 10 “^5 based on previous works 

In our simulations we use conical hoppers with different tilt angles a = 
45°, 60° or 75°, measured from the horizontal plane, and study discharges of 
20000 grains for orifice diameters D = 6R, 8R, 12R, 16R and 24i?. The walls 
are simulated by the force law given by equations and where the overlap 
length ^ is that of the grain with the wall surface. Additionally, i?eff is chosen 
to be equal to the radius of the grain. The other material parameters are equal 
to those of collisions between grains. 

The system is initialized by releasing the grains from a regular cubic lattice 
configuration inside the hopper with the grain centers separated typically by 3R 
or 4:R with small random velocities while the hopper orifice is kept closed. The 
orifice is opened after a time of typically 0.5 s (5 x 10^ time steps) and then the 
discharge begins. 

In order to maintain a stationary flow, grains that are more than one di¬ 
ameter below the outlet are reinserted at randomly chosen positions at the top 
surface of the material in the hopper, this procedure is carried out at regular 
intervals of 0.1 s (10^ time steps). Thus, the number of particles inside the 
hopper is kept approximately constant all along the discharging process. 

A snapshot of the system configuration can be seen in Figure The Carte¬ 
sian and the spherical coordinate system are shown, for both cases the origin 
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is chosen to coincide with the cone vertex. Also, the hopper angle with the 
horizontal plane, a, is defined. 



Figure 1: Snapshot of hopper and flowing grains. Cartesian and spherical coor¬ 
dinates are shown, respectively, as (x^y^z) and 


3 Results and Discussion 

In order to visualize the steady-state flow, we create a velocity field. It is 
obtained by constructing a regular square grid over the plane (p, cp) representing 
by a color in each cell the modulus of the average velocity taken over grains inside 
that cell and the neighboring cells. That is, we average also over 0 and thus, by 
this procedure, we exploit the azimuthal symmetry of the system. Additionally, 
the velocity is averaged over snapshots from the whole simulation, snapshots 
are taken every 0.005s and the duration of each simulation is typically between 
5s and 10s. Some of those fields are shown in Figure for a few choices of 
orifice widths. Note that, by increasing the orifice size, the overall velocity 
increases. The iso-velocity curves, i.e. the borderline between two different 
color shades, are always crossing the entire horizontal extension of the hopper 
for orifice size D = SR, while for D = 16R there are enclosed regions close to 
the wall. Also, the general shape of these curves changes with the orifice size, 
being more horizontal for smaller orifice sizes and increasingly inclined for larger 
ones. This means that for the same distance from the vertex of the hopper the 
average velocity of the grains is larger at the center of the hopper than close to 
its walls. 

The velocity field analysis may be useful to differentiate extreme states such 
as a fully developed funnel flow and a very typical mass flow. However, it 
is unable to quantify more subtle flow differences. Thus, in order to better 
characterize the flow, we contrast the two most different regions of the hopper, 
that is, we compute the velocity profile along p, the radial distance from the 
hopper cone vertex (Figure Q, close to the wall of the hopper, Vy^{p), and along 
its central line, Vc{p). Furthermore, we generalize the concept of Mass Flow 
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Figure 2: Velocity field for a = 60° and D = SR, 16R and 24i?, from left to 
right. 


Index to a Local Mass Flow Index profile, /lmf(p), similarly as in Gonzalez- 
Montellano et al. [^. It can be obtained directly from the velocity profiles 
as 


^lmf(p) = 


^w(p) 

Vcip) ' 


We apply this analysis to data from simulations of discharges in hoppers 
of different orifice sizes, Figure]^ In panel (a), a power law decay is observed 
for the velocity profiles along p. As discussed before, this behavior is not sur¬ 
prising, however, discrepancies between and Vc are visible. We believe that 
these discrepancies carry information about the flow regime and should not be 
disregarded. For small orifice sizes, namely D = 6R and D = SR, the slope 
of the velocity along the center of the hopper appears to be larger than that 
along the walls, until they collapse into approximately the same curve further 
up in the hopper. On the other hand, for large orifice sizes, D = 16R and 
D = 24i?, there is a large gap between Vc and all along p, and the slopes 
are very similar. Therefore, the gap is approximately maintained even for large 
p. One could have analyzed this also through the /lmf in Figure [^. The 
simulations with small opening diameters show /lmf curves increasing with p 
until they reach values close to unity, meaning the wall and center velocities are 
approaching identical values at higher positions. On the other hand, the curves 
for large opening widths fluctuate around much smaller /lmf values, being a 
consequence of the approximately constant gap between and This seems 
to be a precursor of funnel flow where the wall velocity is much smaller than 
the center velocity all along the hopper. 

Additionally, we present the analogous velocity profiles for hoppers of differ¬ 
ent opening angles. Figure [^. The same general behavior is observed, that is, 
a power law for both velocities. Differences between and are larger close 
to the outlet and decrease with increasing p. It is interesting to note that the 
/lmf curves for each angle a, in Figure]^, are identical within the fluctuation 
range. This is a very interesting effect that we intend to investigate further for 
other opening angles. 

According to the literature 


14 23 , the velocity radial profiles should decay 
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Figure 3: Velocity profiles for center (solid line) and wall (dashed line) velocities 
in log-log scale (a) and Local Mass Flow Index (b) for different orifice diameters 
and a = 60°. The profiles are power laws, but slightly different inclinations 
between Uw and Vc can be observed for some cases what causes the I^mf to 
rise. The larger the orifice diameter the larger is the gap between center and 
wall velocities and consequently the lower is /lmf- 


as a power law with exponent equals to 2, as calculated for an incompressible 
fluid in the Introduction of this paper. For all the cases, we obtained exponents 
around that value, the smallest value being 1.868(3) for R = 12 and a = 45° 
and the largest value being 2.176(5) for R = 2A and a = 60°, both for the 
wall velocity profile. However, we observe in Figures and that for the 
center velocity the slope of the profile varies along p. This is easily observed, for 
instance, for R = 12 and a = 45° (lower curves in Figure]^) where the center 
and wall velocity profiles are very close and approximately parallel for large p. 
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Figure 4: Velocity profiles for center (solid line) and wall (dashed line) velocities 
in log-log scale (a) and Local Mass Flow Index (b) for different hopper angles 
and D = 8R. As in figure [^, power law behavior is observed for the velocity 
profiles. Here, the differences between the center and wall velocities close to the 
outlet are larger for smaller 0^ however, the difference vanishes for large p in all 
cases. This is represented in panel (b) by the low values of /lmf for small p 
that then rise up to values close to unity. 


but they separate close to the orifice. We believe this is an effect intrinsically 
related to the flow regime of the discharge and the local density heterogeneities. 
Moreover, the differences between the central and the wall velocity profiles can 
be highlighted by the Local Mass Flow Index. In most cases, the /lmf shows an 
increase with p, presenting intermediate values, around 0.8, close to the outlet 
and rising further up until a value close to unity which means that and Uw 
are equal. This effect is not observable by any other analysis. Also, for the 
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case where the measured /lmf is approximately constant, namely D = 2AR and 
a = 60°, its value is relatively small, that is, close to the limit for funnel flow 
regime, even though it is expected to present mass flow. 

In summary, our results show that the Local Mass Flow Index profile varies 
along p for many hopper angles and opening widths raising concerns over the 
characterization based solely on the Mass Flow Index, where this information 
is lost. The flow regime characterization is far too complex to be reduced to 
a single index. Nevertheless, the analysis of the /lmf profiles was shown to 
be a very powerful tool to describe the more subtle aspects of the flow clearly. 
Although our simulations were limited to mass flow regime, we expect that more 
distinctive features will be observed for discharges in the funnel flow regime, 
what is part of our current investigations, which demands considerable larger 
CPU times. 


4 Conclusions 

We have investigated the influence of the orifice size and hopper opening angle 
on the velocity fields and profiles of spherical grains flowing through a conical 
hopper. 

The fields show the presence of low velocity zones close to the walls of the 
hopper, only observed in systems with large outlet sizes. They seem to be pre¬ 
cursors of stagnation zones that appear in fully developed funnel flow. Also,for 
larger orifice widths as compared to the smaller ones studied, the iso-velocity 
curves present a larger inclination, that is, the average velocity at the center of 
the hopper is larger than the velocity close to its wall for the same distance from 
the hopper vertex. Although these differences do not characterize different flow 
regimes, they show that there is a spectrum of behaviors for each flow regime. 

The velocity profiles along the axis of the hopper and close to the walls were 
studied, as well as their ratio, the /lmf profile. Both velocity profiles show 
approximately a power law dependence with the distance from the vertex of the 
hopper for the considered cases. In hoppers with large orifice widths, there is an 
upshift of the velocity profile along the center of the hopper as compared to that 
close to the walls, but the slopes are roughly the same. This is accompanied 
by a /lmf that fluctuates around a certain value that can be as low as 0.65 for 
the largest orifice width analyzed. On the other hand, for small outlets, the 
difference between the profiles decreases with the distance from the opening. 
This decrease is clearer when the angle of the hopper is changed. The smaller 
the angle of the walls of the hopper with the horizontal plane, the larger the 
slope of the central velocity profile compared to the wall velocity profile. In all 
cases, the velocity profiles do not cross each other, so that, for a distance farther 
from the orifice, the difference of the slopes becomes negligible. This decrease 
in difference between the velocity profiles is seen as an increase of /lmf- The 
/lmf stabilizes with values close to one in the regions where the velocity profiles 
collapse. 

We believe that the presented comparative analysis of the velocity profile at 
the center and close to the walls of the hopper accompanied by the /lmf profile 
is a promising tool for the classification the granular flow regimes. It provides 
a broader and more intelligible picture of the flow state than the techniques 
currently applied in the field and it lacks the set-up sensitivity of other methods. 
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We are currently working on applying this analysis to different silo geometries, 
such as flat and conical bottom cylindrical silos. In those cases a fully developed 
funnel flow is more easily accessible and, for that regime, we expect the /lmf 
analysis to present some new insights. 
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